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Résumé 

Des représentations stochastiques de solutions de problèmes de Dirichlet 
déterministes linéaires et non linéaires sont déduites de l'application de la 
formule de Itô. Ces représentations sont utilisées pour établir des algorithmes 
de calcul par simulation de marches aléatoires. Les méthodes numériques as- 
sociées sont appliquées à des exemples de problèmes linéaires et non linéaires. 
Les résultats des essais avec une fonction source, des estimations de temps 
d'arrêt et des courbes de régression quadratique sont présentés. Des solu- 
tions u du problème de Dirichlet Au = au 3 sont calculées par cette méthode 
purement stochastique pour des valeurs de a positives et négatives. 
Mots clés : Problèmes de Dirichlet déterministes linéaires et non linéaires — 
Équations différentielles stochastiques — Fonctionnelles intégrales — Résolu- 
tion numérique stochastique — Marche aléatoire 



Abstract 

In this paper, we are interested in numerical solution of some linear boundary 
value problems with Dirichlet boundary part, by the means of simulation of 
random walks. We use a probabilistic interprétation of solution u, assum- 
ing that the coefficient and the boundary data are sufficiently smooth, and 
applying Itô's formula. From thèse stochastic représentations of solution, 
we extend some algorithms obtained for standard boundary conditions to 
the quasi-linear source of the type f(u) = au 3 . For positive and négative 
parameter a, we then obtain numerical results by applying the stochastic 
methods based upon thèse generalized algorithms. 

Keywords: Linear and nonlinear Dirichlet BVP; SDE; Probabilistic repré- 
sentation; Numerical stochastic method; Random walks 
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Chapitre 1 
Introduction 



Le calcul des solutions d'une équation différentielle du type u" = f(u(t)) 
est d'un intérêt constant (voir le récent article [3]). Notre objectif est prin- 
cipalement de calculer des solutions de problèmes de Dirichlet u" = au 3 mis 
en exergue dans [16] où a = 1. 

La résolution numérique de problèmes aux limites peut être faite, de 
manière classique, en mettant en œuvre des méthodes de différences finies 
ou d'éléments finis, associées à un maillage de l'espace. 

Des méthodes de Monte-Carlo peuvent également être appliquées ; elles 
conduisent, après discrétisation des équations, à des traitements particuliers 
liés à la géométrie locale des domaines (voir, par exemple, [6], [2], [H] et [9]). 

Cette étude explore une autre voie en utilisant la théorie des processus 
de Markov qui fournit des représentations intégrales pour les solutions des 
problèmes aux limites linéaires et non linéaires déterministes stationnaires. 
On sait, en particulier, que le problème de Dirichlet admet une représentation 
intégrale [5] utilisée pour établir des algorithmes de calcul de la solution par 
simulation de marches aléatoires [Ï5] . 

Les méthodes numériques qui en découlent ne nécessitent pas d'entrer en 
mémoire un maillage de discrétisation dans le cas linéaire. Dans le cas non 
linéaire, on utilise le maillage élémentaire support des marches aléatoires. 
Ceci conduit à une programmation courte et facile à vérifier pas à pas. 

De manière plus précise, on résout un problème de Dirichlet déterministe 
et stationnaire dans des cas linéaires et non linéaires. On montre notamment, 
à partir d'exemples, que l'algorithme stochastique converge plus rapidement 
par relaxation dans le cas non linéaire. 

Le chapitre 2 est consacré aux représentations stochastiques des solu- 
tions de problèmes de Dirichlet linéaires et non linéaires. Les approches des 
représentations par réalisations des processus sont présentées dans le chapitre 
3 et les algorithmes de calcul correspondants sont établis. Des simulations 



numériques des marches aléatoires sont effectuées dans le chapitre 4 pour 
calculer des valeurs des solutions. Une approche fonctionnelle est présentée 
dans le chapitre 5. 



Chapitre 2 

Représentations des solutions 



Ce chapitre est consacré à la représentation des solutions de problèmes 
déterministes de Dirichlet linéaires ou non linéaires par des intégrales sto- 
chastiques. 

Considérons un ouvert borné G de M. d (la dimension d étant fixée). On 
note x G M d la variable d'espace et n la normale unitaire intérieure définie 
sur la frontière dG. 

À partir d'équations différentielles stochastiques (EDS) (pour ces der- 
nières, voir par exemple, [TU [S] et les références associées), nous pouvons 
représenter les solutions de problèmes aux limites en faisant intervenir l'espé- 
rance d'intégrales : fonctionnelles de trajectoires et de fonctions aléatoires 
solutions des EDS. L'application de la formule de Itô en liaison avec les EDS 
conduit à la représentation des solutions. 

Sur le plan fonctionnel, les représentations que nous obtenons n'ont été 
établies en général que sous des hypothèses de régularité sur le domaine G 
et sur les données aux limites (se reporter aux références mentionnées dans 
la suite de ce chapitre). Néanmoins, nous utilisons ici les représentations 
stochastiques lorsque les données ne sont pas nécessairement régulières : la 
frontière dG du domaine peut, par exemple, présenter des coins ou des arêtes 
- cette procédure s'est déjà montrée efficace par exemple dans [SIEJE] pour 
les équations discrétisées et dans [15] pour le problème de Dirichlet. 

Dans le paragraphe 1, nous représentons la solution du problème de Di- 
richlet. Le paragraphe 2 est consacré à des problèmes de Dirichlet quasi- 
linéaires où le second membre dépend de la solution. 



2.1 Problème de Dirichlet 

Considérons le problème de Dirichlet : 

-|Au = / dans G _ _ 

u = g sur <9G 

où u est la fonction inconnue, définie sur G, à valeurs réelles et les données 
sont les fonctions scalaires / et g, définies respectivement sur G et dG. 

Introduisons le processus de Markov : X x = x + W t , t > , où Wt est un 
processus de Wiener standard à valeurs dans M. d tel que l'on ait : X£ = x. 

Soit G le complémentaire de G dans M. d . Le temps d'atteinte de G est 
défini par : 

t x = inf { t > | X x e G } . 

Lorsque la condition : E[r x ] < +oo est vérifiée en tout point x de G, la 
formule de Itô appliquée à u ( X* ) conduit à la représentation de la solution 
de (12. ip . que nous écrivons, pour des raisons techniques, sous la forme : 

u(x) = E[Y], xeG, (2.2) 

avec la variable aléatoire : 

Y= f T f(Xndt + g(X* x ). (2.3) 

Jo 

Des résultats fonctionnels sur ce problème et la représentation associée peuvent, 
par exemple, être trouvés dans [T], [5] et [7]. 

Considérons le temps d'atteinte r et montrons que r dépend du coefficient 
de diffusion a. 

Soit B(0, r) une boule ouverte de M. d , de centre O et de rayon r, de frontière 
C(0,r). 
Considérons u : B(0,r) — > M. la solution unique du problème : 

■^aAu = 1 dans B(0, r) 
u = sur C(0,r) 

où a > est supposé constant. 

Alors u admet la représentation fl2.2H2.3p avec le processus : 



(2.4) 



X? = x + y/âW t . (2.5) 

D'après le raisonnement appliqué par page 253, avec a = 1], on en 
déduit que le temps d'atteinte r x du complémentaire de B(0, r) par le pro- 
cessus X x , à partir d'un point intérieur x, vérifie : 

2 112 

u{x) = E[t x } = r ~ m , xeB(0,r). 

LL (Jj 
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Par conséquent, lorsque a tend vers zéro par valeurs positives, E [r x ] tend 
vers +00. 

De plus, lorsque a décroît, les accroissements du processus (12.51) décrois- 
sent avec a/ô car dXf = y/a dWt ; ceci limite notre champ d'investigation 
lors des essais numériques. C'est pourquoi, avant toute résolution numérique, 
on s'intéresse au comportement de la solution u lorsque a tend vers zéro sur 
G pour éviter les éventuelles difficultés numériques. 

2.2 Problème de Dirichlet quasi- linéaire 

Considérons le problème : 

-7; Au = f(x,u) dans G 
2 V ^ (2.6) 

u = g sur dG 

où / est une fonction donnée de (x, u), f non linéaire en u, supposée régulière. 

Appliquons la méthode des approximations successives : 
A l'ordre m, substituons / par / (u( m ~ r) ) dans (12. ip . alors (I2.2p et (I2.3P 
donnent la représentation de la solution de (12.61) sous la forme itérative : 

u (m) {x) = E[Y^} , xeG, (2.7) 

avec la variable aléatoire : 

Y {m) = f f(x^u( m -^{XÏ)yt + g{X x T ). (2.8) 

Sous certaines conditions, en particulier sur la régularité du second membre 
/, la suite Y^ tend vers la représentation stochastique de la solution u du 
problème 12.61 (voir le chapitre [5]) . 

L'objectif principal de ce travail est le calcul de solutions, supposées suffi- 
samment régulières et bornées, du système non linéaire : 

Au = au 3 dans G =]0; 1[ 
u(0) = (2.9) 

t*(1) = 1 

avec le paramètre a G M et où le second membre / est une fonction non 
linéaire de u. 



2.3 Commentaires 

Les représentations des solutions ont été établies à partir d'équations 
différentielles stochastiques (EDS) (pour ces dernières, voir par exemple, [TU 
[8] et les références associées). 

Sur le plan fonctionnel, les représentations que nous avons obtenues, n'ont 
été établies en général que sous des hypothèses de régularité sur le domaine G 
et sur les données aux limites. Néanmoins, nous utilisons les représentations 
stochastiques lorsque les données ne sont pas nécessairement régulières : la 
frontière dG du domaine peut, par exemple, présenter des coins ou des arêtes 
- cette procédure s'est déjà montrée efficace par exemple dans [SIEJE] pour 
les équations discrétisées et dans [TS] pour le problème de Dirichlet sans 
terme de potentiel (pour cette notion, voir [T3]). 

Les représentations associées aux problèmes de Dirichlet (12.11) peuvent 
être trouvées, par exemple, dans [TJ, [3 page 127] et [TJ pages 244 et 364], où la 
solution est supposée de classe C 2 ( G ). Des résultats fonctionnels concernant 
la régularité de la solution des problèmes de Dirichlet ( 12. ip . et la régularité 
de ses dérivées, sont rappelés dans [2] où les conditions sur les données sont 
telles que la solution est de classe C 2 ( G ). 

Dans le cas non linéaire, des représentations itératives du type (l2.7H2.8p 
peuvent être trouvées dans [5j, mais seulement dans le cas d'une équation 
parabolique semi- linéaire. 



Chapitre 3 

Algorithmes de résolution 



Ce chapitre a pour objet de donner des réalisations des processus associés 
aux représentations du chapitre précédent. Ces réalisations permettent de cal- 
culer effectivement les solutions des problèmes aux limites. Le traitement des 
conditions aux limites se ramène à un comptage des nombres d'absorptions 
sur les frontières du domaine. Le traitement de la non-linéarité conduit à des 
itérations des algorithmes de calcul des solutions de problèmes linéaires. 

3.1 Problème de Dirichlet 

La solution du problème (12. ip est représentée par la moyenne d'une varia- 
ble aléatoire Y fonction d'un processus de Wiener standard [To] . 

La représentation de ce problème par les équations (l2.2H2.3j) . page [7J 
montre qu'il suffit de calculer NT valeurs approchées de Y, notées Y\, . . . , Y NT , 
pour obtenir u(x) : 

NT 

Chaque valeur Y n est obtenue à partir d'une simulation d'un processus de 
Wiener de la manière suivante : le pas h > étant fixé, on simule le processus 
Xf par la suite de vecteurs aléatoires Xç>, ■ ■ ■ , X k , ■ ■ ■ définis par les formules 
de récurrence : 

Xq = x G G , initialisation, 

X k+l = X k + hD k , k G N, 

où D k est un vecteur aléatoire tel que, si e i; i = 1, . . . , d, est la base canonique 
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de M d , alors 

d 

i=l 

vérifie : 

Prob(D fc = e^ = Prob(L> fc = -e*) = — , % = 1, . . . , d . 

2a 

En dimension trois, par exemple, on engendre une suite, nécessairement 
finie, de nombres : 

Uq, . . . , Uk, • • • , Un , 



simulant une variable aléatoire U uniformément distribuée sur l'intervalle 
[0, 1] à l'aide d'un générateur de nombres pseudo-aléatoires ; puis on pose, 
pour tout k : 

l 7 <U k < 1 -^- =► D k = (-1)% , z = 0,...,5, 
oo 

avec 

% 

7 = 1+ - 

J [2 

où les crochets désignent la partie entière de l'argument. 

En dimension deux, à chaque intervalle de temps, le déplacement se fait : 

1. soit à pas constant h dans l'une des deux directions d'un repère or- 
thonormé, direction choisie au hasard avec une probabilité 1/2, et avec 
un des deux sens choisi au hasard avec une probabilité 1/2, comme 
le montrent la figure I3.H -A et l'algorithme 1 du pas élémentaire de la 
figure l372| où la fonction random renvoie un réel appartenant à [0, 1[; 

2. soit à pas constant h\/2 dans l'une des deux bissectrices d'un repère 
orthonormé, direction choisie au hasard avec une probabilité 1/2, et 
avec un des deux sens choisi au hasard avec une probabilité 1/2, comme 
le montrent la figure 13.11 - B et l'algorithme 2 du pas élémentaire de la 
figure 13.21 où la fonction randi renvoie un entier pair ou impair pour 
chaque coordonnée. 



Chaque simulation donne une réalisation X k , < k < N, d'une marche 
aléatoire dans G telle que : 

X.q = x G (j , . . . , X-N—i G Ct et -X-n y 1 " ; 
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(X 1 ,X 2 



H 



H 




Figure 3.1 - Pas élémentaires d'une marche aléatoire dans le plan, avec une 
probabilité de 1/4 chacun 



Algorithme 1 du pas élémentaire 

Ul :=random; U2 : =random ; 

if UK0.5 then 
if U2<0.5 then X2:=X2 + H 
else X2:=X2 - H 
else 

if U2<0.5 then X1:=X1 + H 
else X1:=X1 - H; 



Algorithme 2 utilisé sous MATLAB 



m =[xl;x2] ; 

m =m+h.*(-l) .~randi(2, [2, 1] ) ; 



Figure 3.2 - Algorithmes du pas élémentaire 

dont la trajectoire est la courbe polygonale reliant les points X k et X k+1 , 
k — 0, . . . ,N — 1, successivement. 

Le processus aléatoire X* est alors discrétisé selon la méthode d'Euler ; 
en particulier, si h est le pas d'une marche aléatoire simulée, alors on a les 
accroissements en pas et en temps suivants : 

SX? = h et ôt = h 2 /d. 

À chaque marche aléatoire simulée X k , on définit et on calcule l'encaisse- 
ment Y n : 

N 



Y n = ôt-J2f(Xk)+9i(X N ) 



fc=i 



Puis on répète cette démarche NT fois pour obtenir un échantillon de taille 
NT pour la variable aléatoire Y et on en déduit u(x). 
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L'entier NT correspond au nombre de marches aléatoires simulées, donc 
au nombre d'absorptions par G. Cet entier NT est l'indice maximum de la 
boucle principale du programme. La deuxième et dernière boucle du pro- 
gramme suit pas à pas chaque marche aléatoire en additionnant au fur et à 
mesure la valeur de la fonction source /, puis, pour finir, l'effet de l'absorption 
avant de stocker l'encaissement Y n correspondant. 

Considérons un domaine G de M 2 et un point (X1D,X2D) donné dans G. 
Alors un algorithme possible Al du problème de Dirichlet se présente sous 
la forme : 

YN := 0.0; 
for i:=l to NT do 
begin 
XI := X1D; X2:= X2D; YA:= 0.0; 

while INTER(X1,X2) = true do 
begin 

YA:= YA + F(X1,X2); 
XI := XI + pas élémentaire aléatoire; 
X2:= X2 + pas élémentaire aléatoire; 
end; 

YN:= YN + H * H * YA/2 + G1(X1,X2); 
end; 
U:= YN / NT; 



Figure 3.3 - Boucle du problème de Dirichlet 

La fonction INTER prend la valeur vraie lorsque le point (XI, X2) est dans G 
et fausse sinon. La fonction F définit l'effet de la source en chaque point de la 
marche dans G et Gl définit l'effet d'absorption sur la frontière du domaine. 

3.2 Problème de Dirichlet quasi-linéaire 



La solution du problème (12 .6p est représentée par la moyenne d'une varia- 
ble aléatoire Y fonction d'un processus de Wiener standard. 
La représentation de ce problème par les équations 02.7H2.8p . page [8], montre 
qu'il suffit de calculer NT valeurs approchées de Y, notées Yi, . . . , Y NT , pour 
obtenir u(x) : 



NT 

NT 



"(•'■) TT^S^- 



=i 
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Chaque valeur Y n est obtenue à partir d'une simulation d'une marche aléatoire 
selon l'algorithme 13.21 comme pour la résolution du problème (12. ip . 

Comme indiqué dans le paragraphe I3.1[ chaque simulation donne une 
réalisation X k , < k < N, d'une marche aléatoire dans G, le processus 
aléatoire Xf étant discrétisé selon la méthode d'Euler. 

À chaque marche aléatoire simulée X k , on définit et on calcule l'encaisse- 
ment Y n : 

N 



Y n = St-^2f(u(X k ))+g(X t 



fc=i 

Puis on répète cette démarche NT fois pour obtenir un échantillon de 
taille NT pour la variable aléatoire Y et on en déduit u(x). Dans le cas non 
linéaire, le calcul est effectué aux points de discrétisation du domaine (points 
supports des marches aléatoires), puis réitéré à partir des valeurs obtenues 
en ces points. 

Le cas non linéaire demande quelques explications. Non seulement la solu- 
tion est calculée sur l'ensemble des points de discrétisation du domaine, mais 
encore la non-linéarité impose des itérations des algorithmes précédents. 

Considérons un domaine G = ]0,L[. Alors un algorithme possible A2 du 
problème de Dirichlet non linéaire (12 .6p se présente sous la forme : 

y„ Problème non linéaire avec condition de Dirichlet 
% Subdivision en maxpt intervalles du domaine G = ] 0, L [ 
maxptl = maxpt + 1 ; 

h = L / maxpt ; % Pas de la marche aléatoire 
% miter % Nombre d'itérations 

% nt % Nombre de marches aléatoires par point 

°/ U % Initialisation de U 

% Itération du problème non linéaire 
for iter = limiter 
for J = 2: maxpt 

ys = 0.0 ; dir = 0.0 ; 

for 1=1 : nt % Nb de marches aléatoires par point 
m = J; ya = 0.0; 
while ( Inter(m)==l ) 

ya = ya + F(U(m) ) ; 
m = m + (-l)~randi([0, 1]) ; 
end 

dir = dir + G( m ) ; 
ys = ys + ya ; 
end 
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U(J) = ( h*h*ys + dir ) /nt ; 
end 
end 
°/ Affichage des valeurs obtenues 

La variable maxpt est le nombre d'intervalles de la subdivision de G , 
iter le nombre d'itérations, miter le maximum d'itérations, U le tableau 
d'argument J des valeurs successives de U obtenues à chaque itération. 

La marche aléatoire se fait d'un point x m à un point voisin x m _i ou x m+ i 
de la subdivision du domaine ; le pas aléatoire élémentaire est donc discret : 
on passe d'un indice M à un indice voisin M-l ou M+l, compris entre 1 et 
maxpt + 1. L'entier nt correspond au nombre de marches aléatoires simulées 
à chaque itération, donc au nombre d'absorptions par G à chaque itération. 

En chaque point Xj de la subdivision, une valeur de la solution u(xj) 
est estimée, intégrée à la boucle de subdivision, donc prise en compte pour 
l'estimation des valeurs de la solution aux points suivants x^, k > j. Plus 
précisément, à l'étape m, on calcule une valeur de u™ +l à l'aide des valeurs de 
l'étape m en cours : u™,...,u™, et des valeurs de l'étape précédente m — 1 : 

v m-l m-l 



Dans le paragraphe 14.21 du chapitre suivant, nous nous intéresserons plus 
particulièrement au système (12 .9p et, compte tenu des essais numériques et 
de la convergence numérique, nous serons amenés à modifier cet algorithme. 

3.3 Commentaires 

La vitesse de convergence des algorithmes dépend du pas h de la marche 
aléatoire ; le processus brownien X discrétisé et le temps d'atteinte r de 
la frontière dG en dépendent également. Lors des essais numériques, notre 
champ d'investigation est limité à des valeurs de h qui permettent au pro- 
cessus d'atteindre la frontière en des temps raisonnables, donc aux vitesses 
de calcul des ordinateurs utilisés ; pour éviter d'éventuelles difficultés numé- 
riques, on s'intéresse au comportement des temps d'atteinte avant toute 
résolution numérique. 
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Chapitre 4 
Essais numériques 



Des simulations numériques de la méthode stochastique introduite précé- 
demment sont présentées dans ce chapitre. Ces essais sont destinés à valider 
les représentations du chapitre 2 ainsi que les méthodes approchées du cha- 
pitre 3. Enfin, le système de Dirichlet non linéaire (J2.9P est résolu à l'aide du 
mouvement brownien approchée par des marches aléatoires. 

4.1 Problème de Dirichlet 

Considérons le problème (12.11) dans la couronne G définie par : 
G = Ux,y) G M 2 | 1< yjx 2 + y 2 <?\ 

(cf. figure HZED . avec les données suivantes : 

/ = dans G 

g = 4 sur C(0; 1) = {(x,y) G M 2 | x 2 + y 2 = l} 

g = 6 sur C(0; 3) = {{x,y) G M 2 | x 2 + y 2 = 9}. 

Ce problème de Dirichlet admet pour solution exacte : 

«<*,») = 4 + 2 !îiv^Z. (41) 
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R = 3 




Figure 4.1 - Couronne circulaire G de frontière C(0\ r = 1) U C(0; R = 3) 



À l'aide de l'algorithme Al, on a obtenu, dans un premier temps, les 
résultats du tableau I4.1[ page [T71 avec les données suivantes : 



Point considéré : 
Nombre de tirages 
Valeur exacte : 



(oc, y) = (2;0) 

NT = 10 5 
u(2;0) « 5,26186. 



Table 4.1 - Problème de Dirichlet -- Calcul de la solution en (2; 0) 



Pas d'une marche 


Valeur calculée 


Erreur relative 


h 


u c 


\(u c -u)/u\ 


0,1 


5,27040 


1,6 -HT 3 


0,05 


5,26876 


1,3 -10- 3 


0,025 


5,26730 


1,0 -îo- 3 


0,005 


5,25766 


0,8 -HT 3 



et, dans un second temps, l'évolution de la solution calculée jusqu'à 10 6 
tirages avec h = 10~ 3 de la figure 14.21 où la solution calculée est portée en 
ordonnée tous les 10 4 tirages. 
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Évolution de la solution u en fonction du nombre de tirages 



5, 2653 -r 



5,2619 



5,2596 




2-10 



4 • 10 5 6 • 10 6 

Nombre de tirages 



10 b 



Figure 4.2 - Courbe d'évolution de la solution jusqu'à 10 6 tirages 
théorique : u = 5, 2619 — Pas de la marche simulée : h = 10 -3 



Valeur 



Considérons des couronnes circulaires et rectangulaires de la figure 14.31 de 
la page [T9l 

Des résultats des essais numériques sont représentés dans les figures 14.31 
14.41 et I4T51 ; ils ont été obtenus avec un pas h = 0.1 et un nombre de marches 
aléatoires égal à 5000 par point. 

Dans un premier temps, des évaluations des temps d'atteinte ont été 
effectuées et reportées sous chaque couronne de la figure 14.31 

Dans un second temps, les essais numériques ont conduit aux représen- 
tations de la solution de la figure 14.41 de la page dHJ 

À l'aide de l'outil "Curve Fitting Tool" de MATLAB, l'approximation des 
solutions de la solution conduit aux coefficients suivants (avec un intervalle 
de confiance à 95%) : 

1. Couronne circulaire : f(x) = a + 61n(x)/ln(3), 

a = 4.014 (4.011, 4.016), b = 1.986 (1.981, 1.99) 

2. Couronne rectangulaire 1 : f(x) = pix 2 + P2X + ps, 

pi = -0.4289 (-0.4351, -0.4227), p 2 = 2.531 (2.505,2.556), 
j9 3 = 2.26 (2.237,2.284) 

3. Couronne rectangulaire 2 : f(x) = p\X 2 + p 2 x + p 3 , 

p x = -0.1137 (-0.1201, -0.1072), p 2 = 1.445 (1.418,1.471), 
p3 = 2.673 (2.648,2.698) 
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Circular ring 



Rectangular ring 1 



Rectangular ring 2 






2 


Hitting times 


2 




1.5 


I 


1.5 


1 


^wMtj 


1 


0.5 


/x 


0.5 





Z 7 \J 






1 2 3 

Space variable 




Figure 4.3 - Temps d'atteinte de la frontière de couronne et des parties 
internes et externes 



Dirichlet problems on circular and rectangular rings 




Figure 4.4 - Solutions approchées des trois couronnes 
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Polynomial curves litting of hitting tii 




Figure 4.5 - Temps d'atteinte des frontières des trois couronnes 

Les estimations des temps d'atteinte des frontières des couronnes et leurs 
courbes de régression quadratique sont données sur la figure 14.51 de la page 
I2Ô1 Les courbes de régression admettent pour équations : 

1. Couronne circulaire : f(x) = — 0.9932a; 2 + 3.883a; — 2.76 

2. Couronne rectangulaire 1 : f(x) = — 0.4609a; 2 + 1.6a; — 0.811 

3. Couronne rectangulaire 2 : f(x) = -1.081a; 2 + 4.329a; - 3.244 

Les temps d'atteinte des frontières et des parties intérieures et extérieures 
aux couronnes sont précisés dans l'annexe El 

4.2 Problème de Dirichlet quasi-linéaire 

Dans un premier temps, pour valider l'algorithme A2, résolvons le pro- 
blème (J2.6P dans l'intervalle G = ]0,L[, avec les données suivantes : 



/(«) 



<7i(0) 



u 

2 exp [TT~x 



dans G 



1, ^(L) = (l + L)(l-ln(l + L)) 



où le second membre / est non linéaire en u. 

Ce problème de Dirichlet admet pour solution exacte 



u(x) = (1 + x) (l - ln(l + x)\ . 
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avec les données suivantes 



Longueur de l'intervalle 
Nombre de points : 
Nombre de tirages : 
Pas d'une marche : 
Valeur initiale : 
Nombre d'itérations : 



iltats du tableau 14.21 ont été obtenus 


L 


— 


1 


MAXPT + 1 


= 


51 


NT 


= 


10 3 


h 


= 


2 x io~ 2 


u 


= 





MITER 


= 


10 



(4.2) 



L'intervalle G a été subdivisé en MAXPT = 50 sous-intervalles de même lon- 
gueur. La solution a été calculée, à chaque itération, en chacun des points 
de discrétisation. Après dix itérations, le maximum de l'erreur relative aux 
points de discrétisation a été calculé et vaut : 1,62 x 10~ 2 au point 0,72. 

Table 4.2 - Problème de Dirichlet non linéaire 



Point 


Valeur calculée 


Valeur approchée 


Erreur relative 


x 


u c 


u 


\{u c -u)/u\ 


0,2 


1,0013 


0,9998 


0,15 x 10~ 2 


0,5 


0,8920 


0,8918 


0,02 x 10~ 2 


0,72 


0,79999 


0,78720 


1,62 x 10~ 2 


0,8 


0,7415 


0,7420 


0,07 x 10~ 2 



Dans un second temps, considérons le problème (12.61) . correspondant à 
(J2.9P issu de [16], dans l'intervalle G = ]0, 1[, avec les données suivantes : 

f(u) = — au 3 dans G 

9(0) = 
<7(1) = 1 

Lorsque a = 1, ce problème de Dirichlet admet une solution qui vérifie [Ï6] : 

Vx G G, x( 1+v ^) /2 < u{x) < x . (4.3) 

À l'aide de l'algorithme AI (mais sans relaxation), les résultats de la 
figure 14.61 ont été obtenus avec les données suivantes : 
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Nombre de points : 


MAXPT + 1 = 


= 51 


Nombre de tirages : 


NT -- 


= 2 x 10 3 


Pas d'une marche : 


h = 


= 2 x kt 2 


Valeur initiale : 


u = 


E 



Courbe représentative de la solution calculée u c à l'itération 39 



u. 



1 




y(x)=x 



}y ,yU=x {l+ ^ )/2 



1 



U c 



G = ]0,1[ 

Figure 4.6 - Représentation de la solution du problème de Dirichlet non 
linéaire : u" = u 3 , u(0) = 0, u(l) = 1 en trait continu 

Le nombre d'itérations MITER a été remplacé par le test d'arrêt : 

x\ < u{xi) < Xi est réalisé (cf. (14. 3p ) et m(xj_i) < u(xi), pour tout i = 

2,. . ., MAX, c'est-à-dire que la solution calculée aux points Xi de discrétisation 
est bien encadrée et est strictement croissante. Sur la figure 14.61 le test est 
satisfait à l'itération 39. 

Avec le logiciel MATLAB, le programme prend la forme suivante : 

°/ Problème non linéaire avec condition de Dirichlet 

clear ail 

rng( 'shuff le' , 'v5unif orm') ; 

tic 

% Domaine G = ] 0, 1 [ 

1 = 1.0 ; 

% Subdivision en maxpt intervalles 

maxpt = 20 ; 
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maxptl = maxpt + 1 ; 

h = 1 / maxpt ; 7» Pas de la marche aléatoire 
°/ Nombre d'itérations 
miter = 500 ; 

°/ Nombre de marches aléatoires 
nt = 20000 ; 
% Initialisation de U 
U = : h : 1 ; 

°/ Fonctions encadrantes et lissage 
x = : maxpt ; 

y = (h.*x).~(0.5*(l+sqrt(5))); 
z = h.* x ; 
hspline = 0.25*h; 

% Itération du problème non linéaire 
for iter = limiter 
for J = 2: maxpt 

ys = 0.0 ; dir = 0.0 ; 

for I = 1 : nt % Nb de marches aléatoires par point 
m = J; ya = 0.0; 
while ( (m > 1) && (m < maxptl ) ) 

ya = ya - 0.5* U( m )* U( m )* U( m ) ; 
m = m + (-l)~randi([0, 1]) ; 
end 

if m > 1 dir = dir + 1.0 ; end 
ys = ys + ya ; 
end 

U(J) = ( h*h*ys + dir ) /nt ; 
end 

figure (iter) hold ail cla 
xx = : hspline: maxpt ; 
yy = spline(x,U,xx) ; 
plot (xx,yy) 
plot(x,y) ;plot(x,z) 
end 
toc 

Pour la valeur a = 1, nous avons obtenu une courbe représentative de la solu- 
tion calculée numériquement par la méthode stochastique et l'erreur relative 
par rapport à des calculs par différences finies suivant la fonction MAT- 
LAB 'bvpc' appliquée à de nombreuses valeurs initiales dans la commande 
MATLAB 'bvpinit'. Des courbes sont données dans les figures ETT1 où l'erreur 
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Nonlinear Dirichlet Pi 
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r □( Slochastic solution 


versus Interpolaled 
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\ / 


V--^ 


N/ 




-0.005 










- 


-0.0, 


\ 








- 


-0.015 


-V 








- 



2/5 1/2 3/5 7/10 4/5 

Space variable 



Figure 4.7 - Solution pour a = 1 - Erreur relative empirique 



relative empirique atteint un maximum inférieur à 2% au voisinage de zéro 
et où la méthode stochastique est régularisante par rapport à la méthode aux 
différences finies (à condition de prendre des nombres de marches aléatoires 
et d'itérations suffisamment grands, ce qui a un coût en temps de calcul). 

Compte tenu des difficultés de convergence numérique et du temps de 
calcul, le programme de calcul a été amélioré en premier lieu par le choix 
aléatoire du point de départ des marches aléatoires, ce qui évite le biais 
du choix croissant des points, et par l'introduction d'un calcul de moyennes 
progressives à l'aide du poids bary qui permet de cumuler l'effet des marches 
aléatoires successives. Par rapport au programme précédent, le programme ci- 
après a donné une convergence numérique plus rapide, en particulier pour a = 
— 1 (voir la figure fl~8l où l'évolution de l'erreur relative empirique (u — x)/x 
est reportée). 

% Itération du problème non linéaire 
miterl = miter * maxpt ; 
for iter = l:miterl 

m = randi( [ 2 , maxpt ] ) ; % Choix aléatoire 
ys = 0.0 ; dir = 0.0 ; 

for 1 = 1 : nt °/ Nb de marches aléatoires par point 
indice = m ; ya = 0.0; 
while ( (indice > 1) && (indice < maxpt 1 ) ) 

ya = ya + U3demi( indice ) ; % + au lieu de - 
indice = indice + (-1) ~randi( [0, 1]) ; 
end 

dir = dir + (indice>l) ; 
ys = ys + ya ; 
end 

U(m) = ( (h2*ys + dir )/nt + bary(m)*U(m) )/(bary(m)+l) ; 
U3demi(m) =0.5 .* U(m).* U(m).* U(m) ; 
bary (m) = bary (m) + 1 ; 
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onlinear Diric 


let Problem u"= - u 3 








An = -«■" 






















Figure 4.8 - Solution pour a = — 1 - Courbes de convergence numérique 



end 



Ce programme a donné les résultats représentés dans les figures 14.81 

Pour faciliter l'écriture d'un programme de calcul parallèle avec MATLAB, 

nous effectuons le changement de variable suivant : 



ku 



avec k = sgn(a) 



Alors le système devient : 

Av 
v{0) 

v(l) 



sgn (a) v 



sgn (a) J 



dans G 



(4.4) 



Soit M > assez grand. On cherche une solution v bornée par M, dérivable à 
l'ordre deux, à dérivées continues et bornées sur l'intervalle fermé G* = [0, 1]. 
Alors on a : 

veCl(G*)^f(v)eCl(G*) 

où f(v) = sgn(a) v 3 . 

La fonction / est lipschitzienne en v, uniformément sur G, ce qui satisfait 
les conditions du Théorème 15.11 du chapitre [5j 

Ces propriétés sont présentées dans l'annexe \M 

La méthode a conduit aux résultats de la figure 1431 avec une subdivision 
en 20 intervalles et un nombre de marches aléatoires égal à 20 000 en moyenne 
par point. 

Les valeurs de a correspondant à u'(l) = Oet àw'(l) = —1 ont été 
obtenues avec la méthode stochastique et les fonctions bvp4c de Matlab et 
NDSolve de M athematica (voir le tableau 14. 3|) . 

Pour ces deux cas, la méthode stochastique a conduit aux résultats de la 
figure 14.101 avec une subdivision en 20 intervalles et un nombre de marches 
aléatoires égal à 10 000 en moyenne par point. 
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Nonlinear Dirichlet problems: u" = a u 




1/10 1/5 3/10 2/5 1/2 3/5 7/10 
Space variable 



Figure 4.9 - Solutions pour a e -4, -3, -2, -1, 1, 2, 3, 4, 5 
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Table 4.3 - Valeurs de a pour lesquelles u'{\) = et u'(l) 



u'(l) 


Méthode 
stochastique 


MATLAB 
bvp4c 


MATHEMATICA 
NDSolve 



-1 


a = -3.43755656 
a = -4.39 


a = -3.43761588 
a = -4.33564372 


a = -3.43759287 
a = -4.33553126 



0.6 



0.2 



Nonlinear Dirichlet problems: u" = a u 



1/10 



1/5 



3/10 



2/5 1/2 3/5 

Space variable 



7/10 











J< /^ 




//^<^^ 














Jp^ I I I I I I I 


I I 



4/5 



9/10 



Figure 4.10 - Solutions pour u'(l) = et u'(ï) = -1 

4.3 Commentaires 

Pour a > 0, on trouve des solutions très proches (moins de 1%) de celles 
obtenues avec les logiciels Mathematica (commande NDSolve, en particulier 
avec les méthodes "Shooting" et " ExplicitRungeKutta" ) et MATLAB (en 
particulier avec la commande bvp5c). Les courbes ont donc la même forme 
convexe que celle de la figure IA.1I de l'annexe IA1 

Pour a < 0, on trouve des solutions concaves pour a > —4.5. Par contre, 
a = —4.5 et a = —5 conduisent à une convergence numérique "lente" dans un 
premier temps (par accumulation de marches aléatoires), puis à une diver- 
gence numérique de manière "explosive" dans un second temps (où la solution 
atteinte numériquement est prise comme valeur initiale). Par conséquent, il 
ne faut pas se limiter à une convergence numérique, mais il est nécessaire 
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de réitérer la méthode en prenant la valeur obtenue dans un premier temps 
comme valeur initiale; le "bruit blanc" du mouvement brownien permet de 
tester la stabilité de cette valeur obtenue lors de la première étape ; si cette 
valeur est stable, alors une solution est atteinte ; sinon, la méthode ne permet 
pas d'affirmer que la valeur obtenue dans un premier temps soit une solution 
(même instable). 

Pour a = —3.43755656 et a = —4.39, les courbes obtenues ont les formes 
des figures E]2] de l'annexe iBl La forme IÂ731 de l'annexe iBl n'a jamais été ob- 
tenue lors des essais numériques de la méthode stochastique, sauf de manière 
transitoire. La méthode stochastique privilégie les solutions les "plus stables" 
en un sens qui reste à définir. Les autres solutions sont soit instables, soit 
métastables quand elles possèdent un domaine d'attraction réduit, si ce do- 
maine existe (les calculs n'ont pas permis de trouver de tels domaines). 
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Chapitre 5 

Approche fonctionnelle 



5.1 Introduction 

Le problème de Dirichlet non linéaire (I2.6P a été résolu numériquement à 
l'aide de la représentation f l2.7H2.8p par reports successifs des valeurs obte- 
nues. Ce chapitre présente une approche fonctionnelle de ce type de représentation 
sous des hypothèses supplémentaires d'existence et de régularité de la solu- 
tion. Cette approche par une méthode de type Picard est utilisée dans [5], 
exclusivement à l'aide de la formule de Feynman-Kac, pour établir certains 
résultats concernant la solution de problèmes de Cauchy dans M. d avec des 
conditions particulières sur les données et sur les opérateurs. L'objectif prin- 
cipal de ce chapitre est de justifier, autant que possible, l'utilisation des 
représentations qui permettent de résoudre numériquement le problème non 
linéaire (12. 6p . 

Nous considérons le domaine G = ]0,L[ où L > est une constante 
donnée ; nous notons x la variable d'espace appartenant à G ; n désigne la 
normale unitaire intérieure définie sur dG. 

Considérons le problème : 

— ~ Au = f(u) dans G 

(5.1) 
u = g sur dG 

où l'inconnue u est définie de G dans M et les données sont : 

/ : R — >R 

u i — > f(u) 
g : dG — y R 

x i — y g(x) 
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Une représentation stochastique permet l'application d'un théorème de 
point fixe. Ceci conduit à montrer dans le paragraphe 15.2^ l'existence et 
l'unicité de la solution u de (15 .ip avec un second membre non linéaire f(u) 
sous des conditions sur /. 

5.2 Une méthode de point fixe 

Sur le plan fonctionnel, il s'agit de se ramener à traiter des représentations 
stochastiques de la solution de (15. ip par une méthode de point fixe. 

La suite de ce paragraphe concerne principalement la résolution du sys- 
tème (15. ip avec / suffisamment régulière, par utilisation d'un théorème de 
point fixe de type Picard. 

On introduit les espaces fonctionnels suivants : 

(i) C£ (G), respectivement C£ ( dG), l'espace des fonctions bornées, k 
fois dérivables en i 6 G, respectivement x G dG, à dérivées bornées ; 

(U) Cb(G), l'espace des fonctions continues et bornées, muni de la norme 
sup : 

u H- \u\oo = max{ \u(x)\ : xEG}. (5.2) 

Les étapes de la méthode peuvent être définies ainsi : 

1. on prend «W dans C%(G) et on vérifie que l'on a : f(u^) dans 
C 2 b (G); 

2. en supposant que, pour / ( «W ) dans C\> ( G ), il existe une solution u^ 
du problème : 

-iAîif 2 ' = /(««) dans G 

,, (5-3) 

m (2) = g sur dG 

vérifiant : U® eC b (G), 

on obtient une représentation stochastique associée ; 

3. on définit l'application u^ = Q [vS 1 '} et on montre que G a un point 
fixe unique dans un espace adapté. 

On doit vérifier préalablement que, pour tout u G Cl ( G ), on a / {u) G 
Cl {G), pour pouvoir définir la représentation stochastique associée au pro- 
blème. 
Exemple. 

La fonction f{u) = au 3 est de classe C°° en u et vérifie : 

ueC 2 b (G)^f(u)eC 2 b (G). (5.4) 
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La seconde étape nécessite la mise en évidence d'une formulation stochas- 
tique associée à ( 15.3p . noté plus simplement : 



(5.5) 



— 7; Au = f dans G 
u = g sur dG 

où les fonctions f et g vérifient, par hypothèse, les propriétés : 

/ e C 2 b (G) (5.6) 

9 G C 2 b {dG) (5.7) 

Proposition : On suppose /, g donnés vérifiant (j5.6H5.7ft . 

Alors, le système linéaire (15.51) admet une solution unique u dans C\ ( G). 
Démonstration : voir [10J. 

Il reste à établir le résultat de cette approche fonctionnelle. Nous adaptons 
au cadre stochastique, à l'aide des hypothèses d'existence et de régularité, 
les arguments utilisés dans le cadre déterministe. 

Théorème 5.1. On suppose que : 

1. le domaine en espace est défini par : 

G = ]0,L[, avec < L < +00 ; 

2. le second membre f vérifie : 

et la fonction f est lipschitzienne en u, uniformément sur G ; 

3. la fonction g vérifie \5.7\) ; 

4- les données du système linéaire 115. 5\) sont telles que, pour tout f dans 



Cb(G), il existe une solution unique u telle que : 

ueC h {G). (5.8) 

Alors A5.1\) admet une solution unique u dans Cb(G), qui est limite de fonc- 
tions dans Cl (G) . 

Cette solution u est représentée, pour tout x G G , par : 



u(x) = E 



f(u(x:),x: 



+ E[g{X-)\ 

(5.9) 
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Démonstration : Pour obtenir une représentation stochastique de la so- 
lution, nous restons dans un cadre continu et borné. 

Soit vS 1 ' dans C& (G ). Considérons la fonction / définie sur G par : / = 
/ (m*- 1 ) ). Les hypothèses de la proposition 15.41 entraînent : 

/ = /(««) eCa(G), 

puis, d'après l'hypothèse (15.8p . on obtient l'existence et l'unicité de la solu- 
tion, notée «(% de (15. 3p où l'on a posé : / = / (u^ ). 

Lorsque l'on prend u^ 1 ' dans C% (G), alors d'après la proposition 15.2^ on 
a : u^ G C 2 ( G ), et, en appliquant la formule de Itô avec une démonstration 
similaire à celle du théorème 5.1 de [5j pages 167-168], on obtient la représen- 
tation suivante : 



u {2 \x) 
De (15.8p . l'application : 



E J T f(uU(X*),Xl 



+ E[g(XÏ)] (5.10) 



g : C b (G) — ► C b (G) 

(i) __>. „( 2 >=0( W (D) 



tr ' i — >■ u 



où u^ 2 ^ est la solution unique de (I5.3P avec / = / ( u^ ) , est bien définie. 

Montrons que Q est une contraction sur Ct,(G). Soit u^\ vP-> dans 
Cb(G) et les solutions correspondantes dans Cf, ( G ) : u' 2 ) = <? ( m 1 - 1 ^ ) , m 1 - 2 ^ = 

Nous utilisons la notation générique : ôq = q — q . En particulier, notons : 

ÔU W = u (l) _^(l) ) 5m (2) = u (2) _^(2) ; 5/ ( M (l)) = /( W W) -/(««), 

5uM = sg ( u w ) = s («w ) - e (# ). 

Cherchons un réel ^ dans ]0, 1[ tel que : 

\sg(u {1) )\ <u\ôu {1) \ . 

I \ / I oo — I I oo 

Les solutions u^ et û^ de (I5.3P sont telles que leur différence ôu^ vérifie : 

-|A(W 2 )) = (6f)(uU) (G) 
Ôu&> = (dG) . 

D'après l'hypothèse (15.8p . ce système admet une solution unique : 

ôu (2) eC b (G) . 
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Alors, pour tout a; G G, on a : 

(ôu {2) )(x) = E 



{5f)(u^{X:),Xl 



IJQ 



Cette représentation est donnée, par exemple, dans [5] pour une solution 
dans Cl {G) d'un problème linéaire. Comme on a seulement ôu^ G C& ( G), 
cette représentation est déduite de la formule de Itô, comme dans [2], en 
approchant Su^ par une suite de fonctions régulières dans Cf (G). 
La représentation précédente entraîne pour tout x G G : 



6u®i 



X) 



< E 



\(ôf){ u w(x:),x:)\ds 



Dans la suite de cette démonstration, M désigne une constante générique 
telle que : < M < +oo . 

Comme la fonction / est lipschitzienne en u, uniformément sur G, on a : 



\ôu {2) (x)\ < ME 
On en déduit : 

\ôu {2) (x)\ < ME 
et : 



(<5m (1) )(X s *)| ds 



\ôu {1) \ ds 



WxeG. 



VxeG. 



5u {2 \x) I < ME[t] \5u {1) \ 



VxeG 



où E[t] < +oo. 

Pour terminer la preuve de l'existence d'un réel v dans ]0, 1[ rendant l'ap- 
plication Q contractante sur Cb(G), nous raisonnons par récurrence. 

A l'ordre m = 2, on prend l'inégalité précédente. 



A l'ordre m, on fait l'hypothèse de récurrence suivante : 

\6u {1) \ , WxeG. 



faW(I) i < [Mm 

1 (771—1)! 



A l'ordre m + 1, on a successivement : 
1. de la même façon qu'à l'ordre m = 2 

\ôu {m+1 \x)\ < ME 



(ôu {m) )(X*)\ ds 



33 



2. d'après l'hypothèse de récurrence à l'ordre m : 

\m—l 



\ôu {m+1 \x)\ < ME 
3. par calcul de l'intégrale : 



ME[t] 



5u {m+l \x) 
ôu {m+1) (x) 
ôu {m+1) (x) 



< 



< 



< 



M r ' 



(m-1)! 

M m 
(m-1)! 

(ME[t] 



(m-1) 

\ôu w \ E 
\ôu {1) \ E\t 

I I oo L 

m 

- \5uW\ 



\5u (1) \ ds 



E[t 



m—l 



ds 



VxeG. 



ml 



La récurrence étant établie, on en déduit : 

{ME[r\y 



8u {m+1) 



< 



SuM 



ml 



(5.1i; 



Or la série \_. 



MEW] 



m=l 



ml 



converge, donc il résulte de f)5.1ip que (u^ m ' ) m>1 



forme une suite de Cauchy qui converge dans l'espace Cb(G) vers une limite 
unique u. 

Par récurrence, nous définissons l'itérée de Q à l'ordre m, notée Ç( m \ Pour 
m suffisamment grand, G^ est contractante de coefficient v tel que : 

I V / 1 00 — I 1 00 

L'application contractante Q^ 1 admet donc un point fixe unique dans Cj, ( G). 
D'après un corollaire du théorème du point fixe de Banach, on sait que Q 
admet aussi un point fixe unique u dans Cb(G) tel que : u = Ç{u) . 

Finalement, on en déduit que u est la solution unique de (15. ip . et que u 
est limite de fonctions de Cl ( G ), représentées sous la forme (I5.10p . 

La représentation (15.91) résulte de (I5.1UJ) et de u = G(u) . □ 

5.3 Commentaires 

Les résultats généraux d'Analyse Fonctionnelle appliqués dans ce chapitre 
peuvent être trouvés, par exemple, dans [TU! H]. 

Le résultat du théorème 15.11 peut s'étendre à des conditions aux limites 
de Neumann ou de Robin avec un ouvert borné G de M. d , dont la frontière dG 
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et la normale intérieure n sont suffisamment régulières (voir les conditions 
nécessaires à l'existence de ce genre de représentations stochastiques dans [SJ 
pages 166 et 254] et [H page 117] où dG, respectivement n, est de classe C 3 , 
resp. C 2 ). 
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Chapitre 6 
Conclusion 



Des méthodes stochastiques de calcul de solutions de problèmes détermi- 
nistes de Dirichlet linéaires et non linéaires ont été présentées. Les représenta- 
tions résultent d'une application rigoureuse de la formule de Itô pour les 
problèmes réguliers et de son application formelle pour les problèmes à do- 
maines non réguliers. Les cas linéaires et quasi-linéaires peuvent ainsi être 
traités. Ces méthodes sont basées sur des représentations stochastiques qui 
donnent directement des algorithmes aisément programmables. 

Les programmes sont courts, faciles à construire et à vérifier pas à pas. 
On évite d'entrer en mémoire un maillage de discrétisation du domaine et de 
gérer les tableaux de numérotation qui l'accompagnent. 

Suivant le problème considéré (géométrie du domaine, expressions de la 
fonction source / et des conditions aux limites g), le temps de calcul est plus 
ou moins long avec le processeur séquentiel que nous utilisons. 

Comme les méthodes de Monte-Carlo classiques, cette méthode stochas- 
tique admet une vitesse de convergence en 1/y/NT. Les problèmes quasi- 
linéaires nécessitent davantage de tirages. Des propriétés essentielles des 
méthodes de Monte-Carlo sont conservées : d'une part, le calcul de la so- 
lution en un point choisi indépendamment des autres dans le cas linéaire, 
d'autre part, l'adaptation au calcul parallèle dans les cas linéaires et non- 
linéaires. Mais notons également que les marches aléatoires sont simulées 
simplement à partir d'épreuves répétées de type Bernouilli. 

La résolution numérique basée sur des représentations stochastiques peut 
s'appliquer en dimension deux ou trois sans difficultés autres que celles liées 
au temps de calcul. Les représentations stochastiques peuvent aussi être as- 
sociées à la méthode itérative de Picard (voir le chapitre [5] pour des résultats 
partiels d'existence et d'unicité de la solution). 

Sur le plan fonctionnel, ces méthodes stochastiques ont nécessité des pro- 
priétés supplémentaires de régularité, en particulier la dérivabilité de la solu- 
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tion au sens classique, et des propriétés sur les données telles que le système 
linéaire ( 15. 5 p admette une solution dans Cb(G). Concernant la régularité 
des solutions de problèmes aux limites non linéaires ou de problèmes posés 
dans des ouverts à frontière non régulière, signalons que dans pH] et [5] pour 
des problèmes paraboliques non linéaires et linéaires, des estimations et des 
représentations stochastiques sont établies pour des solutions suffisamment 
régulières, de classe C 2 en espace. 

Remerciements Ce travail a bénéficié du soutien de l'Université de La 
Réunion. 
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Annexe A 

Système non linéaire Au = au , 

iz(0) = 0, u(ï) = 1 



Considérons le système non linéaire (12. 9p : 

An = au 3 dans G =]0; 1[ 

u(0) = (A.l) 

u(l) = 1 

avec ael. 

Lorsque a = 0, ce problème admet la solution évidente u(x) = x. 

Lorsque a = 1, ce problème de Dirichlet admet une solution qui vérifie [16] : 



Vx G G, X ( 1+VE V 2 < u{x) < 



X. 



A.l Changement de variable et propriétés 

Le changement de variable : 

v = ku avec k = sgn(a)y|a| 

donne le système suivant : 

Av = sgn(a) v 3 dans G 

v{0) = (A.2) 

v(l) = sgn(a)y|a| 



3S 



On suppose que v est bornée par M > assez grand. On s'intéresse à une 
solution v bornée par M, dérivable à l'ordre deux, à dérivées continues et 
bornées sur [0,1]. 

Proposition A.l. SoitC^ l'ensemble des fonctions définies sur [0, 1], dérivables 
à l'ordre deux, à dérivées continues et bornées. Pour tout a G M, on a : 

veC 2 b ([0A})^f(v)eCU[0,l}) 

où f(v) = sgn(a) v 3 . 

Preuve : évidente. 

Pour satisfaire les conditions du Théorème 15.11 du chapitre El montrons 
la proposition suivante. 

Proposition A. 2. La fonction f définie par f(v) = sgn(a)v 3 est lipschit- 
zienne en v, uniformément sur G. 

Preuve : Soient v±, Vi G C 2 ( [0, 1] ). 

I f(vi) - f(v 2 ) | < |sign(a)| | v? - vf | 

< max I v\ + v it> 2 + v ? I x I Vi ~ Vo I 

~ [0,1] ' ' 



< 3 M 2 I vi - v 2 



a 



Proposition A. 3. Soit v une fonction continue bornée solution de AA.2\) . 

On suppose que v n'est pas identiquement nulle sur un intervalle [0, A], 
< A < 1. 

1. f" est continue et bornée sur [0, 1]. 

2. Si a > 0, v' ne s'annule pas sur [0,1] et v est strictement croissante 

sur [0, 1] 

3. Si a < et s'il existe x e]0, 1] tel que v'(x ) = 0, alors on a : \v(xq)\ = 



V2\v'(x Q )\. 
Preuve : L'équation v" = sgn(a) v 3 entraîne que v" appartient à Cf 



b ■ 



Intégrons f" = sgn(a) v sur [0, 1]. On a : 



v'v" = sgn(a) v'v 3 
\{v'f = \sgn(a)v' + C 
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où C est une constante réelle. 

Les conditions aux limites entraînent : 

v(0) = =» C = ^'(0)) 2 
u(l) = sgn(a)v1âi =^ C = -(f'(l)) 2 - ^sgn(a) a 2 
Par conséquent, on obtient : 

(vf = ~sgn(a) «* + MO)) 2 (A.3) 

S'il existe x G]0, 1] tel que u'^o) = 0, alors on a : 

^sgn(aMx ) 4 + K(0)) 2 = (A.4) 

Si a > et v non identiquement nulle à droite de zéro, alors il n'existe 
pas Xq annulant v', qui conserve ainsi le même signe sur [0, 1] ; or v(0) = et 
v(l) = 1, donc v' > sur [0, 1] et v est strictement croissante sur [0,1]. De 
plus, v est convexe sur [0, 1]. 

Si a < et s'il existe x e]0, 1] tel que v'(xq) = 0, alors x vérifie : 



v{x f = 2t/(0) 2 
v(x ) 2 = V2\v'(Q)\ 



\v{x )\ = yJV2\v'{0)\ 



D 



Remarque A.l. Comportement de la solution au voisinage de zéro. 
1. Le calcul approché de 

// x i v(h) -v(0) , v(/i) 
v'(0) = lim -^- — ^ = lim -V^ 

fc-H) h h-+Q h 



permet d'estimer : v(xq) = ±y y/2\v'(0)\. 
2. Quelques conséquences de IjA.J^ : 
(a) (v>(l)) 2 -(v>(0)) 2 =hgn(a)a* 



1„2 



(b) Sia>0, alors \v'\ > \v'(0)\ sur ]0, 1] et v'{lf > \a 2 

(c) Si a < et s'il existe x e]0, 1] tel que v'(x ) = 0, alors \v(x )\ > 1 



si et seulement si |i>'(0)| > ^-. 



2 
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(d) Si a < 0, alors \v'\ < \v'(0)\ sur ]0,1], v'(l) 2 = v'(0) 2 - \a 2 et 
v'(l) 2 < v'(0) 2 . De plus, on a : a 2 = 2 (V(0) 2 - v'(l) 2 ). Le calcul 
approché de a tel que v'(l) = 0, permet de calculer une valeur 
approchée de v'(0) 2 , soit : 

t/(l) = =» v'(Of = £ . 

3. Formule de Taylor de v en zéro. 

h 5 
v(h) = hv'(0) + —sign{a)v'(Qf + o(h 5 ). 

ZA} 



A. 2 Formes des courbes représentatives 



Les solutions u du problème (lA.lj) obtenues par la méthode stochastique 
ont été comparées aux solutions obtenues par différences finies. Pour a < 0, la 
méthode stochastique fournit et privilégie une solution strictement positive. 

Indiquons les formes des courbes représentatives de solutions u du pro- 
blème ( lA.ip suivant le signe de a. Pour a < 0, nous ne présentons que les 
formes de courbes où la solution ne s'annule au plus qu'une fois. 

Cas : a > 

u est positive, croissante, continue, bornée et convexe sur [0, 1]. 




Figure A.l - a = 10 

Cas : a < 

Les différentes figures "régulières" obtenues avec le logiciel MATLAB ont 
la forme des figures IA31 et IA.31 où le segment représente la fonction identité, 
pour les cas suivants : 

1. Forme Fi : a — —3 avec 

bvpinit(linspace(0,l,10) , [ 0.5 0.5 ]) 

u est positive, strictement croissante, continue, bornée et concave sur 
[0, 1]. u atteint son maximum en 1. u, u' et u" ne s'annulent pas sur 

]0,1[. 
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Figure A. 2 - Formes F\ et F 2 obtenues avec a < 




Figure A. 3 - Forme F 3 obtenue avec a < 

2. Forme F 2 : a = —3 avec 

bvpinit(linspace(0,l,10) , [ 3 3 ]) 

u est positive, continue, bornée et concave sur [0,1], croissante, puis 
décroissante, u admet un maximum sur [0, 1]. m et u n ne s'annulent pas 
sur ]0, 1[. u' s'annule une seule fois sur ]0, 1[. 

3. Forme F 3 : a = —3 avec 

bvpinit(linspace(0,l,10) , [ -3-3 ]) 

u est continue, bornée et convexe sur [0, 1], décroissante, puis croissante. 
u admet un minimum sur [0, 1] et s'annule une seule fois sur ]0, 1[. u" 
s'annule sur ]0, 1[ en même temps que u, qui est d'abord convexe, puis 
concave, u' s'annule une seule fois sur ]0, 1[. 

Formes de courbes obtenues à l'aide de la commande bvpc5 de MATLAB 
pour a = — 1 selon les choix de valeurs a priori de la solution et de sa dérivée 
en dix points : 

1. Les valeurs bvpinit (linspace(0, 1, 10) , [ 1 1 ]) donnent la même 
courbe que celle obtenue avec la méthode stochastique, à moins de 0.3% 
près. Cette courbe est la plus proche de celle de la courbe identité. 
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Obtained curves for a=-1 




0.4 0.5 0.6 

Space variable 



Figure A. 4 - Courbes obtenues avec a 



•1 



2. Les valeurs bvpinit(linspace(0, 1, 10) , [ 1.79 1 ]) donnent une 
courbe convexe avec un minimum proche de —4 atteint en x ~ 0.45. 

3. Les valeurs bvpinit (linspace (0, 1, 10) , [ 1.8 1 ]) donnent une 
courbe oscillante de type sinusoïdal où la fonction s'annule 2 fois. Le 
maximum atteint en x ~ 0.275 vaut environ 7.6 et le minimum atteint 
en x ~ 0.75 vaut —7.6 environ. La courbe devient polygonale aux 
voisinages des extrema. 

4. Les valeurs bvpinit (linspace(0, 1, 10) , [ 1.801 1 ]) donnent une 
courbe convexe identique à celle obtenue avec 
bvpinit(linspace(0,l,10) , [ 1.79 1 ]). 

5. Les valeurs bvpinit (linspace(0, 1, 10) , [ 1.81 1 ]) donnent une 
courbe polygonale (irrégulière) où la fonction s'annule 7 fois. L'algo- 
rithme de la méthode des différences finies ne converge pas dans ce 
cas. Cette courbe à oscillations de très grandes amplitudes n'est pas 
représentée sur la figure IA.4I 

6. Les valeurs bvpinit (linspace(0, 1, 10) , [ 2 1 ]) donnent une cour- 
be concave. Le minimum atteint en x ~ 0.625 vaut environ 3.25. 

Remarque A. 2. La méthode stochastique ne fournit que la première courbe; 
numériquement, les valeurs des autres solutions calculées par une méthode 
de différences finies ont été prises comme valeurs initiales de la méthode 
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stochastique, ce qui a conduit soit à une divergence (dernier cas), soit à une 
convergence vers la première courbe (autre cas) ; en un sens qui est à préciser, 
la méthode stochastique donne une solution stable (les effets de "bruit blanc" 
ne font pas diverger l'algorithme pour certaines valeurs a priorij. 

Pour a = —4, les quatre courbes suivantes ont été obtenues à l'aide de 
MATLAB : 



Obtained curves for a=-4 




20 
Space variable 



Figure A. 5 - Courbes obtenues avec a 



-4 



avec successivement : 

bvpinit(linspace(0,l,10) , [ 0.5 1 ]), 
bvpinit(linspace(0,l,10) , [ 11]), 
bvpinit(linspace(0,l,10) , [ -1 -1 ] ) et 
bvpinit(linspace(0,l,10) , [ -1.5 -1 ]). 

Dans ce dernier cas, la solution oscillante de type sinusoïdal est différente 
de celle obtenue avec a = — 1. 

Seule la solution du premier cas a été obtenue avec la méthode stochas- 
tique. 

Pour a < —5, aucune solution n'est obtenue à l'aide de la méthode sto- 
chastique. Pour a = —5 et a = —6, les solutions calculées avec MATLAB 
ont été utilisées comme valeurs initiales de l'algorithme stochastique et ont 
conduit à la divergence numérique de la méthode stochastique (qui donnent 
des résultats de plus en plus grands jusqu'à la limite du logiciel). 
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Remarque A. 3. Pour une même valeur de paramètre a < 0, les différentes 
solutions calculées avec MATLAB dépendent des valeurs a priori que nous 
avons choisies tandis que V éventuelle solution calculée avec la méthode sto- 
chastique est unique et ne dépend pas des valeurs initiales que nous avons 
utilisées. 
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Annexe B 

Temps d'atteinte et leurs 
courbes de régression 

B.l Introduction 

Etant donnée une partie 1^ de R", le temps d'atteinte de y à partir du 
point x par le processus x + W est la variable aléatoire réelle t v définie par : 

r v : Q x — ► R + 

u i— ► T v (uj) = ini{s>0\X s (uj) eV} 

où Q x est l'ensemble des trajectoires issues de x. 

Par la suite, le temps d'atteinte de V est la fonction Ty : 



Ty = E [Ty] : G 



x i — ► E [t } 



Lorsqu'il n'y a pas d'ambiguité, la partie V n'est pas précisée et Ty s'écrit r. 
D'une façon générale, considérons le temps d'atteinte r et montrons que 
r dépend du coefficient de diffusion D. 

Soit B(0, r) une boule ouverte de M. d , de centre O et de rayon r, de frontière 
C(0,r). 
Considérons u : B(0,r) — y M la solution unique du problème : 



(B.l) 



-\D&u = 1 dans 5(0, r) 
u = sur C(0,r) 

où D > est supposé constant. 

Alors u admet la représentation (j2.2H2.3p avec le processus : 

Xf = x + VDW t . (B.2) 
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D'après le raisonnement appliqué par [7] page 253, avec D = 1], on en 
déduit que le temps d'atteinte r x du complémentaire de B(0,r) par le pro- 
cessus Xf , à partir d'un point intérieur x, vérifie : 

2 I 12 

u{x) = E[r x ] = r Zp , ieB(0,r). 

Par conséquent, lorsque D tend vers zéro par valeurs positives, E 1 [r 11 ] tend 
vers +oo. 

De plus, lorsque D décroît, les accroissements du processus (1B.2J) décrois- 



sent avec \/D car dXf = y D dW t ; ceci limite notre champ d'investigation 
lors des essais numériques. C'est pourquoi, avant toute résolution numérique, 
on s'intéresse au comportement de la solution u lorsque D tend vers zéro sur 
G pour éviter les éventuelles difficultés numériques. 
L'étude du système (12. 9p écrit sous la forme : 

±±Aw = u 3 dansG=]0;l[ 
u(0) = (B.3) 

u(l) = 1 



conduit à considérer un coefficient de diffusion D — - qui tend vers zéro 
lorsque a tend vers +oo. À la limite, la solution du système est la fonc- 
tion indicatrice xi du point 1 ; la solution du système (12. 9p pour a ^ 1 est 
régularisante de xi- 

La vitesse de convergence des algorithmes dépend du pas h de la marche 
aléatoire ; le processus brownien X discrétisé et le temps d'atteinte r de 
la frontière dG en dépendent également. Lors des essais numériques, notre 
champ d'investigation est limité à des valeurs de h qui permettent au pro- 
cessus d'atteindre la frontière en des temps raisonnables, donc aux vitesses de 
calcul des ordinateurs utilisés ; pour éviter d'éventuelles difficultés numériques, 
on s'intéresse au comportement des temps d'atteinte avant toute résolution 
numérique. 

B.2 Résultats des essais numériques : temps 
d'atteinte et courbes de régression 

Les courbes de régression quadratique sont déduites des équations données 
par l'outil Curve Fitting Tool du logiciel MATLAB. 
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Hitting times 




9/5 2 11/5 12/5 13/5 14/5 
Space variable 



Figure B.l - Couronne circulaire - Temps d'atteinte 

B.2.1 Temps d'atteinte r de la frontière de la couronne 
circulaire 

On s'intéresse au temps d'atteinte restreint à un segment porté par un 
rayon, en particulier le temps d'atteinte restreint à l'intervalle [1, 3] sur l'axe 
des abscisses. 

Equation de régression quadratique du temps d'atteinte r de la frontière 
d(0,l)UC 2 (0,3): 

f( x ) = -0.9932x 2 + 3.883a; - 2.76 

Les coefficients sont donnés avec un intervalle de confiance à 95% : 

pi = -0.9932 (-1.007,-0.9795) 
p2 = 3.883 (3.828,3.938) 

p3 = -2.76 (-2.812,-2.707) 

Le temps moyen en seconde pour atteindre la frontière est E = 0.7052 où E 
est la moyenne des E[r Xi ], (xi), i — 1,. . .,p, est une subdivision du segment 
[0, 1] et r Xi est le temps d'atteinte d'un mouvement brownien issu de %{. 
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Equation de régression quadratique du temps d'atteinte du disque interne 
fermé D 1 = V(0,l)UC 1 (0,1) : 

f( x ) = -0.3473a; 2 + 2.064a; - 1.706 

Les coefficients sont donnés avec un intervalle de confiance à 95% : 

pi = -0.3473 (-0.3599,-0.3347) 
p2 = 2.064 (2.013,2.115) 

P 3 = -1.706 (-1.754,-1.658) 

Le pourcentage de marches aléatoires atteignant le disque interne Di est : 
/i =40.8521%= 0.408521. 

Le temps moyen en seconde pour atteindre le disque interne par les 
marches aléatoires absorbées par le disque interne (le temps n'est pas pondéré 
par le nombre de fois où les marches aléatoires atteignent le disque interne 
Di) : Ei = 0.9169 où E\ est la moyenne des -E'prf 1 ], rf l est le temps d'atteinte 
d'un mouvement brownien issu de Xi et atteignant le disque interne Di. 

La pondération donne le temps moyen : ï\ = - Y^=i ^[ T i i ] xn i = 0.6491 
où rii est le nombre de fois où la marche aléatoire atteint le disque interne 
Di à partir du point :tj. 

Equation de régression quadratique du temps d'atteinte du domaine com- 
plémentaire R 2 \ V(0, 3) de V(0, 3) : 

f( x ) = -0.3437a; 2 + 0.6926a; + 1.022 

Les coefficients sont donnés avec un intervalle de confiance à 95% : 

pi = -0.3437 (-0.3539,-0.3336) 
p2 = 0.6926 (0.6517,0.7334) 

p3 = 1.022 (0.9831,1.06) 

Le temps moyen en seconde pour atteindre la partie externe au disque de 
rayon 3 par les marches aléatoires absorbées par cette partie (le temps n'est 
pas pondéré par le nombre de fois où les marches aléatoires atteignent la 
partie externe) : E 2 = 0.9210 où E 2 est la moyenne des E\r 2 % t%* est le 
temps d'atteinte d'un mouvement brownien issu de Xi et atteignant la partie 
externe au disque de rayon 3. 

La pondération donne le temps moyen : T 2 — - Y^=i ^[ T 2 i ] xn i = 0.7438 
où rii est le nombre de fois où la marche aléatoire atteint la partie externe à 
partir du point x t . 

Compte tenu du pourcentage de marches aléatoires atteignant le disque 
interne Di, on retrouve le temps moyen d'atteinte de la frontière : 

E = j\Ti + (1 - fi)T 2 = 0.7052 
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Hitting times 




1 6/5 7/5 8/5 9/5 2 11/5 12/5 13/5 14/5 3 

Space variable 



Figure B.2 - Couronne rectangulaire 1 - Temps d'atteinte 

B.2.2 Temps d'atteinte de la frontière de la couronne 
rectangulaire 1 

Considérons la couronne rectangulaire : G = G\ \ G2 où G\ et G2 sont 
les rectangles : 



G 1 



{ (x,y) G 



max 



(1x1,12/1) < 3} 



G 2 = {(x,î/)Gl 2 max(\x\,\y\)< l}. 

On s'intéresse au temps d'atteinte restreint à l'intervalle [1, 3] porté par l'axe 
des abscisses. 

Equation de régression quadratique du temps d'atteinte de la frontière de 
G : 

f(x) = -0.4609x 2 + 1.6s - 0.811 

Les coefficients sont donnés avec un intervalle de confiance à 95% : 

pi = -0.4609 (-0.487,-0.4349) 
p2 = 1.6 (1.495,1.705) 

P 3 = -0.811 (-0.9106,-0.7114) 

Le temps moyen en seconde pour atteindre la frontière est E = 0.3892. 
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Equation de régression quadratique du temps d'atteinte du rectangle 
fermé [— 1, +1] 2 : 

/(s) = -0.2805a; 2 + 1.535s - 1.175 

Les coefficients sont donnés avec un intervalle de confiance à 95% : 

pi = -0.2805 (-0.2978,-0.2632) 
p2 = 1.535 (1.466,1.605) 

P 3 = -1.706 (-1.241,-1.109) 

Le pourcentage de marches aléatoires atteignant le rectangle interne est : 
/i = 27.5465%= 0.275465. 

Le temps moyen en seconde pour atteindre le rectangle interne par les 
marches aléatoires absorbées (le temps n'est pas pondéré par le nombre de 
fois où les marches aléatoires atteignent le rectangle interne) : E\ = 0.6648 
où Ei est la moyenne des -E'frf*], rf l est le temps d'atteinte d'un mouvement 
brownien issu de Sj et atteignant le rectangle interne. 

La pondération donne le temps moyen : ï\ = - Y^h=i -^t 7 ^] x n i = 0.4447 
où p est le nombre de points de la subdivision (s«), ni est le nombre de fois 
où la marche aléatoire atteint le disque interne D 1 à partir du point Sj. 

Equation de régression quadratique du temps d'atteinte du domaine com- 
plémentaire M 2 \] - 3, +3[ 2 de d : 

f( x ) = -0.1469x 2 + 0.1429s + 0.7954 

Les coefficients sont donnés avec un intervalle de confiance à 95% : 

pi = -0.1469 (-0.1542,-0.1396) 
p2 = 0.1429 (0.1135,0.1724) 

p3 = 0.7954 (0.7675,0.8233) 

Le temps moyen en seconde pour atteindre la partie externe par les marches 
aléatoires absorbées par cette partie (le temps n'est pas pondéré par le 
nombre de fois où les marches aléatoires atteignent la partie externe) : E 2 = 
0.4488 où E 2 est la moyenne des ^[rf*], rf*, rf* est le temps d'atteinte d'un 
mouvement brownien issu de s, et atteignant la partie externe. 

La pondération donne le temps moyen : T 2 = - X]f=i -^t 7 "?] x n i = 0.3681 
où p le nombre de points Sj d'une subdivision du segment [0, 1], n« est le 
nombre de fois où la marche aléatoire atteint la partie externe à partir du 
point Sj. 

Compte tenu du pourcentage de marches aléatoires atteignant le rectangle 
interne [— 1, +1] 2 , on retrouve le temps moyen d'atteinte de la frontière : 

E = fiTi + (1 - fi)T 2 = 0.3892 
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Figure B.3 - Couronne rectangulaire 2 - Temps d'atteinte 

B.2.3 Temps d'atteinte de la frontière de la couronne 
rectangulaire 2 

Considérons la couronne rectangulaire : G = G\ \ Gi où G\ et G2 sont 
les rectangles : 

G, = {(x,y)eR 2 \ sup(|x|,M)<3} 
G 2 = {(x,y)ER 2 \ sup(|x|,M)<l}. 

On s'intéresse au temps d'atteinte restreint à l'intervalle [1, 3] porté par l'axe 
des abscisses. 

Equation de régression quadratique du temps d'atteinte de la frontière de 
G : 



m 



-1.081x 2 + 4.329x- 3.244 



Les coefficients sont donnés avec un intervalle de confiance à 95% : 

pi = -1.081 (-1.087,-1.075) 
P 2 = 4.329 (4.304,4.354) 

P 3 = -3.244 (-3.267,-3.22) 

Le temps moyen en seconde pour atteindre la frontière est E = 0.7324. 
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Equation de régression quadratique du temps d'atteinte du rectangle 
fermé G2 '■ 

f(x) = -0.2849x 2 + 1.86a; - 1.58 

Les coefficients sont donnés avec un intervalle de confiance à 95% : 

pi = -0.2849 (-0.2961,-0.2738) 
p2 = 1.86 (1.815, 1.905) 

p3 = -1.58 (-1.623,-1.537) 

Le pourcentage de marches aléatoires atteignant le rectangle interne G2 
est : /1 = 46.5136%= 0.465136. 

Le temps moyen en seconde pour atteindre le rectangle interne G2 par 
les marches aléatoires absorbées par G 2 (le temps n'est pas pondéré par le 
nombre de fois où les marches aléatoires atteignent le rectangle interne G2) '■ 
Ei = 0.9026 où Ei = Y7i=i -^bi^L P ^ e nombre de points Xi d'une subdivision 
du segment [0, 1], rf* est le temps d'atteinte d'un mouvement brownien issu 
de Xi et atteignant le rectangle interne G^- 

La pondération donne le temps moyen : T\ = Y^ï=i ^[ T i l ] x n i = 0.6501 
où rii est le nombre de fois où la marche aléatoire atteint le disque interne 
G 2 à partir du point Xi. 

Equation de régression quadratique du temps d'atteinte du domaine com- 
plémentaire Il a; Il 1 > 3 de Gi : 

f( x ) = -0.3865a; 2 + 0.7979a; + 1.106 

Les coefficients sont donnés avec un intervalle de confiance à 95% : 

pi = -0.3865 (-0.3981,-0.375) 
p2 = 0.7979 (0.7514,0.8445) 

p3 = 1.106 (1.062,1.15) 

Le temps moyen en seconde pour atteindre la partie externe au rectangle 
G\ par les marches aléatoires absorbées par cette partie (le temps n'est pas 
pondéré par le nombre de fois où les marches aléatoires atteignent la partie 
externe) : E 2 = 1.0275 où E 2 est la moyenne des -E^t^] et rf' est le temps 
d'atteinte d'un mouvement brownien issu de Xi et atteignant la partie externe 
au rectangle G\. 

La pondération donne le temps moyen : T 2 — - 5^f=i ^[ T 2 i ] x n i = 0.8040 
où p est le nombre de points Xi d'une subdivision du segment [0, 1] et rii est 
le nombre de fois où la marche aléatoire atteint la partie externe à partir du 
point Xi. 

Compte tenu du pourcentage de marches aléatoires atteignant le rectangle 
interne G2, on retrouve le temps moyen d'atteinte de la frontière : 

E = f\Ti + (1 - fi)T 2 = 0.7324 
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